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ABSTRACT 

We present a revised method for simultaneous determination of the pattern 
speed (f2p) and star formation timescale (tsr) of spiral galaxies, which is orig- 
inally proposed in our previous work. As this method utilizes offsets between 
molecular and young-stellar arms, we refer to it as the "Offset Method" . Details 
of the method, its application, and results for CO and H« images of 13 nearby 
spiral galaxies are described here. CO data are from our observations with the 
Nobeyama Millimeter Array for 2 galaxies, and from the BIMA SONG for the 
rest. Out of 13 galaxies, we were able to derive Qp and tsF for 5 galaxies. We 
categorize them as "C" galaxies as their offsets are clear. Our findings from 
these galaxies are as follows. (1) The corotation radius calculated by the derived 
Qp is close to the edge of the CO data, and is about half of the optical radius 
for 3 galaxies. (2) The derived tsF is roughly consistent with the free-fall time 
of typical molecular clouds, which indicates that the gravitational instability is 
the dominant mechanism triggering star formation in spiral arms. (3) The tsF is 
found to be almost independent of surface density of molecular gas, metallicity, 
or spiral arm strengths. The number of "C" galaxies and the quality of CO 
data, however, are not enough to confirm these relationships. We also find that 
2 other galaxies show no offsets between CO and Ha, although their arms are 
clearly traced, and categorize them as "N" galaxies. The presence of a bar could 
account for this feature, since these 2 galaxies are both barred. With one galaxy 
excluded from our analysis due to its poor rotation curve, offsets of the remaining 
5 galaxies are found to be ambiguous. Either their dependence on the rotational 
frequency cannot be explained by our picture, or the number or quality of data 
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is not sufficient for tlie analysis. We categorize tfiem as "A" galaxies. The pos- 
sible reasons for this ambiguity are (1) the density wave is weaker, and/or (2) 
observational resolution and sensitivity are not enough to detect the spiral arms 
and their offsets clearly. The former is supported by our finding that the arm 
strengths of "A" galaxies are slightly weaker than that of "C" galaxies. 

Subject headings: galaxies: fundamental parameters — galaxies: individual (NGC 
0628, NGC 3184, NGC 3938, NGC 4254, NGC 4303, NGC 4321, NGC 4535, NGC 
4736, NGC 5194, NGC 5248, NGC 5457, NGC 6181, NGC 6946) — galaxies: spi- 
ral — ISM: Hll regions — ISM: molecules 



1. Introduction 
1.1. Pattern Speed 



Ever since the spiral density wave theory was proposed by iLin fc Shul (119641 ) to solve 
the winding problem of spiral arms, it has been known as the most successful theory to 
explain observational features of spiral galaxies. The pattern speed (Op), defined as the 
angular rotational velocity of a spiral pattern or underlying gravitational potential, is one 
of fundamental parameters in the density wave theory, since it determines the existence and 
location of kinematical resonances. 

The one that has been extensively discussed is the "Corotation Resonance" (CR), de- 
scribed as 

^{r = Rcr) = ^p, (1) 

where Q is the angular rotational velocity of materials. This equation means that at the 
CR, the rotational speed of materials and the pattern are the same. The radius -Rcr, which 
satisfies this equation, is called the "Corotation Radius" . As the rot ation speed of gas in the 



patte rn-rest frame becomes small around the CR, a galactic shock (jFujimotdll968l : [Roberts 



19691 ) would not occur. Since galactic-scale star formation is thought to be enhanced by 
this shock, star formation should be less efficient around the CR. As the relative velocity 
of materials viewed from the pattern changes its sign at the CR, the direction of streaming 
motions due to the galactic shock should be different between the inside and outside of the 
CR. These observable features have been used to locate the CR as described later in this 
section. 

Another important resonance is the "Lindblad Resonance" (LR), which is defined by 

nir = /?lr) = ^^p ± -, (2) 
m 
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where m is any integer but a frequency of the strongest modes or the number of arms is 
mostly used, and k is the epicychc frequency, which is expressed as 



2n\ 1 



(3) 



r dn 

With an assumption that spiral arms are tightly winding, the dispersion r elation for the 
stellar disk is different from that of a gaseous disk (iBinney &: Tremaindll987l ). so that gases 
and stars are theoretically predicted to show different behaviors at the resonances. 

As described above, the pattern speed not only locates the resonances, but also influences 
the star formation activities. Therefore, its determination is very important to the study 
of galaxies. It cannot be determined directly from observations, however, since the pattern 
structure is not a material structure but a density wave. 

Several methods have been proposed for its determination. The earlier and common 
technique is t o locat e specific resonances on radii where some properties of arms change. 
Roberts et al.l (119751 ) adopted flp for 24 galaxies so that the corotation radius is nearly 
coincident with the radius where arni s of H ll regions disappear, or the radial extent of the 
"easily visible" disk. lElmegreen et al.l (119921 ) used 5-band images to locate 5 resonances for 
18 galaxies. Such techniques are, however, subject to uncertainties regarding the sensitivity 
and/or resolution of the imaging observations. In addition, the tracers used must be carefully 
considered, since stars and gases can behave differently at resonances and their distributions 
are not identical. 



Cepa fc Beckmaru (119901 ) derived the star formation efficiency (SFE) in the arm and 
interarm region for NGC 628 and NGC 3992 from Ha and Hi data. The arm-to-interarm 
ratio of the SFE is expected to be larger than 1 as the star formation is thought to be 
enhanced in arm regions, but they found that its value drops to almost unity at a certain 
radius. Since star formation at spiral arms are thought to be less efficient around the CR, 
they concluded that this radius is the corotation radius. It is difficult to derive the SFE values 
in interarm regions, however, since there are fewer stars and gases than in arm regions, and 
thus sensitive observations are required. 

For our Galaxy, a time perio dicity of the star formation history of about 0.5 Gyr has been 

found at the solar neighborhood (IHernandez et al.ll2000l : Ide la Fuente Marcos fc de la Fuente Marcos 
20041 ). Assuming that this periodicity is due to the pa ssage of a spira l potential and that 
the potential has a two-armed pattern as suggested by iDrimmell (120001 ) . the pattern speed 
can be calculated as Qp = fl{r = Rq) — 7r/(0.5 Gyr) = 21 km s~^ kpc^^ A consistent 
result from nume rical simulatio i is of t he stellar and gaseous response to the spiral potential 
are presented by iMartos et al.l (120041 ). These results for the Galaxy also indicate a strong 
relationship between the star formation and the spiral potential. 
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Canzianl (Il993l ) showed that the residual velocity field, obtained by subtracting the ax- 
isymmetric component from the observed velocity field, should be different inside and outside 
the corotation radius, and thus a change in direc tion of the noii - axisy mme tric component of 
the ve locity field could be used to locate the CR. ISempere et al.l ( 119951 ) and lCanzian fc Allen 
( I1997I ) applied this method to Hi and Ha data of NGC 4321, respectively, and they obtained 
a comparable location of the CR. Since this method requires precise kinematical information 
over the entire disk, a velocity field covering the whole galaxy with high velocity and spatial 
resolution is essential. 



Tremaine &: WeinbergI (jl984j ) presented a method which did not use any morphological 
locations of resonances, but used the continuity equation for a surface brightness of galaxies. 
This method is called the "Tremaine- Weinberg method" (hereafter TW method). As it is 
based on the continuity equation, it has been a pplied mostly to early type g alaxies (e.g.. 



Merrifield k KuiikenI 1 19951 : ICerssen et all Il999 



Debattista fc Williams! l2004j ). whose star 



formation is less active than that of late type galaxies. In recent studies, however. Hi and 
CO data have come to be used for the application to gas- domina t ed, la te type galaxies (e.g., 
Westpfahi nJ : Izimmer et all 120041 ). Meanwhile, IPebattistal f |2008l ) showed by A^-body 
simulations that the results of the TW method is s ensitiv e to the uncertainty in adopted 
position angle (P.A.) of outer disk. iRand fc WalhnI (120041 ) applied the TW method to CO 
data of nearby galaxies with varying P.A.s and confirmed that the uncertainty in P.A. made 
the error of the resultant f2p larger. They also showed that the best-fitted f2p only to the inner 
bar was larger than th at of the reniaining disk, even though the galaxy in interest was not 
classified as SB galaxy. iMeidt et al.l (120081 ) used a modified method which explicitly allowed 
a radial variation of fip, and found multiple spiral patterns for the nearby spiral galaxy M 51 
(NGC 5194) with outwardly decreasing Vt-p. The dependence on P.A.s, however, was found 
to be still large. 

The results of nume rical simulations a lso show a large dependence of Jlp on kinematics 
and spiral structure (e.g.. lWada et al.lll998[). so that they c an be used to derive the best value 



of f2p by comparison with observations. lOey et al.l (120031 ) adopted an evolutional model of 
the Hll luminosity function to draw isochrones of massive stars, and fitted these isochrones 
to an observed distribution of Hll regions to derive fip. These approaches, of course, depend 
on their modelings, and often have difficulty estimating the accuracy of the derived values 
and the effect of other parameters on the results. 
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1.2. Star Formation Timescale 



Since all tlie ptiysical processes from molecular clouds to young stars are not yet clearly 
understood, especially for massive objects, it has been difficult to derive a timescale for star 
formation from a theoretical point of view. A number of numerical simu lations regarding the 



Vazquez-Semadeni et al. 



GMC formation a nd successive star formation have been performed (e.g., 
20071 : IWadall2008l ) and have inferred on the timescale, but the full processes are not traced 
since the dynamic range of the size scale and the number of processes are too large to be in- 
cluded in calculations. Instead, it is common to estimate ages of Hll regions from equivalent 
widths of emission lines from ionized gas, such as Ha and H/3, on the basis of population 
synthesis models, which do not require any information about the parental molecular clouds. 



Leitherer et al.l (119991 ) compiled physical processes and parameters for star formation 



such as the initial mass function, star formation history, and metallicity, to provide the 
package Starburst99 to calculate observab le properties such as colors, spectra, equivalent 
widths, and luminosity. iBastian et al.l (120051 ) estimated the equivalent width of Ha emissions 
and derived ages of star cluster complexes in NGC 5194 as 5-8 Myr using this package assum- 
ing solar metallicity. Even if twice the solar metallicity is a ssumed, the derived ages do not 



chang e largely. They also applied the GALEV SSP models (lAnders fc Fritze-v. Alvensleben 



20031 ) to optical colors of individual star clusters and found the ages of the r najority of the 



cluste rs to be 4-10 Myr. With spectroscopic observations of the Br7 line, lGrosb0l et al. 



( 120061 ) derived ages of i^-band knots, presumably sites for massive star formation, to be 
about 7-10 Myr for the nearby spiral galaxy, NGC 2997. They also found that these knots 
were located slightly upstream from the smooth peak of i^-band flux, regarded as the po- 
tential minimum of the spiral arm. 

We should note that ages derived from models depend on parameters used in the model, 
and that the exact age is often difficult to determine due to degeneracies between internal 
extinction, mass, and metallicity. In addition, this age only represents the timescale from 
a newly formed star, and thus has no information about the molecular cloud collapse in 
contrast with the timescale we take into account in this paper. 

As described above, it is in principle difficult to determine the pattern speed and star 
formation timescale from observations, although they are both important parameters for 
understanding kinematics and star formation activities in spiral galaxies. In this paper, we 
propose a method to determine both parameters simultaneously. We refer to this method as 
the "Offset Method" as it uses offsets between CO and Ha arms. This method is a revised 



version of what we have presented in a previous paper (lEgusa et al.ll2004l ). To apply this 
method, we have obtained CO, Ha, and rotation curve data from our observations and the 
literature. 
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We describe the basic idea and characteristics of the method in ^ and processes of the 
apphcation in ^ followed by results in §4l which includes a categorization of the sample, 
discussions for individual galaxies, and properties of each category. We discuss derived 
parameters and the validity of the offset method in §21 and conclude this paper by giving a 
summary in ^ 



2. Offset Method 



This "Offset Method" is proposed for simultaneous determination of the pattern speed 
and star formation timescale by the use of offsets observed in spiral galaxies. A s the basics 
of this method have been presented by our previous paper (lEgusa et al.l l2004j ) . we repeat 
important concepts and explain several modifications here. 



2.1. Basic Idea and Formulation 



For simplicity, we assume that the spiral pattern is rigid, and that materials in the disk 
rotate in pure circular orbits. The former assumption means that the spiral structure is sus- 
tained by the density wave, and that the pattern speed (fip) is constant. The velocity fields 
of most disk galaxies show a spider diagram pattern, which is a velocity field that results from 
pure circular rotations. The streaming motion and the velocity dispers ion generates some 



non-circular motions, but they are typically about 10 -50 km s ^ (e.g., lAdler &: Westpfahl 



19961 : ICombes fc BecquaertI 119971 : iKuno fc Nakail 119971 ) , which is small compared to the cir- 
cular rotational velocity of around 200 km s~^. These observational results support the 
latter assumption. We should note that we do not consider nor include any bar structures in 
our analysis, sinc e the pattern speed of the bar would be different from that of spiral (e.g., 
Wada et al.lll998l ) and particles within the bar potential move in elliptical orbits with higher 
eccentricity. 

We define the star formation timescale, tsF, as the average time for the massive star 
formation from molecular clouds in spiral arms, which are agglomerated or compressed by 
the spiral structure. If the physical processes of star formation at spiral arms do not vary 
extremely with radius, this timescale can be regarded as a constant parameter representing 
a typical value of the entire disk. We should note again that this timescale tsF would be 
different from ages of Hll regions mentioned in §1.2[ since tsF includes the time needed for 
clouds to be collapsed. 



Figure [T] illustrates our concept regarding the relationship between the spiral structure 
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and star formation processes inside the CR, where f2 > fip. If we observe a face-on spiral 
galaxy at t = (the left panel), at t = tsF the same galaxy will be observed as in the right 
panel and the offset distance between the arm of young stars and molecular clouds, d, can 
be written as 

where v is the velocity of materials, and vp is the velocity of the pattern. If the outside of the 
CR is taken into account, massive stellar arms will be seen on the concave side of molecular 
arms, and d will be negative. 

Dividing both sides of equation (jl]) by radius r [kpc], we obtain 



9 



Lp 



— ) [km kpc-^], (5) 



9 = 0.586 



-kms~^kpc~^/ Vkm s ^ kpc ^, 

where Q = v/r, Qp =vp/r, and 9 is the azimuthal offset. This equation shows the relation 
between two observables, Q and 9, and we can rewrite it as 

( ^ ] - ( ]] X ( ] [degree]. (6) 

Vkm s~^ kpc^^/ Vkm s~^ kpc~^ / J VlO^ yr/ 

Assuming tsF to be constant over the spiral disk, ^ is a linear function of Q, since Qp is 
also assumed to be constant. Therefore, by plotting 9 against Q and fitting them with a 
line, both Qp and tsF can be simultaneously determined as a horizontal-axis intercept and 
gradient of the fitted line, respectively. 

The fitting is performed by the least method and each data point is weighted by 
errors in both axes, and AQ. The inclusion of the uncertainty in Q and the use of the 
least method for the fitting are major methodological revisions from what we presented 
in the previous paper. 



2.2. Uniqueness and Merits 

One major advantage of the offset method is its independency from the tight-winding 
approximation for spiral structure and any models for cloud collapse and population synthe- 
sis. The results thus will provide constraints or parameters to theories and models of both 
galactic dynamics and star formation on the basis of observational data. 

Another uniqueness is that it gives Qp and tsF simultaneously, as described in the pre- 
vious section. Measuring offsets in a wide range of radii makes this possible. In addition, 
statistical uncertainties of both Qp and tsF are determinable by the fitting. We should 
emphasize that this is the first method proposed that derives tsF observationally with quan- 
titative estimates of statistical errors. 
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2.3. Requirements for Application 

Application of the offset method, requires tracing both molecular and young-stellar 
arms at a wide range of radii. Nearby galaxies with small inclination angle satisfy this 
requirement, while distant and highly-inclined galaxies do not. Grand-design galaxies are 
usually preferable to see spiral structures. We do not pay attention to whether a galaxy has 
a bar in its center or not, since we put importance on spiral structures only and exclude bar 
regions from the analysis, if any. Images used in the analysis should have sufficient quality 
to detect and resolve spiral arms. For our apphcation, CO and Ha images are used to trace 
molecular gas and young stars, respectively. 

In addition to the image maps, we need information about the rotational velocity. As 
the line-of-sight velocities of an almost face-on galaxy do not give velocities parallel to the 
disk, the inclination angle must not be close to zero. The velocity data are available from 
spectroscopic observations in the optical and radio wavelengths. 

In short, we need nearby spiral galaxies which are moderately inclined and show clear 
CO and Ha spiral structures. To satisfy these requirements, we observed two nearby spiral 
galaxies, NGC 4254 and NGC 6181, in the ^200(1-0) hue with the Nobeyam a Millimeter 



Array (NMA). For NGC 4254, we have combined the single dish data from iKuno et al. 



( 120071 ) with the NMA data to recover the missing flux. Details of the observations and data 



reduction will be pre sented in another p aper. In addition to our data, we selected 11 galaxies 



from BIMA SONG ( iHelfer et al.l l2003l ) . The whole sample of 13 nearby spiral galaxies is 
described in ^3.1[ 



3. Data Analysis 

3.1. Property of Sample Galaxies 

In order to apply the offset method, we selected 13 nearby spiral galaxies. The selection 
criteria are that the galaxy shows a molecular spiral structure in CO data, and that its 
inclination angle is small but not zero. As we pay attention to spiral structures of molecular 
gas and Hl l regions in the analysis, we do not take into account the bar type nor arm class 



(AC) from lElmegreen &: Elmegreenl (119871 ) for the selection. For these sample galaxies, we 
obtained Ha images and rotation curve data from the literature. Most of the images were 
obtained via the NASA/IPAC Extragalactic Database (NED). 

The general properties of our sample are hsted in Table [H Out of the 13 galaxies, 6 are 
classifled as SA (no bar), another 6 are classifled as SAB (weakly barred), and 1 is classifled 
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as SB (strongly barred). Regarding the AC, 12 are grand-design galaxies (AC=9 or 12), and 
1 is a flocculent galaxy (NGC 4736, AC=3). 

We show CO contours overlaid on an Ha image for the sample galaxies in Figure [2J 
Properties of the data are listed in Table [2l The spatial resolution of CO data is about 
3"- 7", and typically 500 pc at the galaxy's distance. As it is comparable to the typical value 
of offsets and arm widths, the CO data quality should be sufficient to resolve offsets for most 
of the sample galaxies. The Ha image seeing is typically about 2", which is smaller than 
that of CO. 



3.2. Rotation Curve 



In general, rotation curves (RCs) are obtained from spectroscopy, such as optical slit 
observations and radio mapping observations. With the CO data cube, we made velocity 
field maps and applied the task GAL in the AIPS software in order to derive the RCs. 

This task fits the tilted-ring model to the velocity field, we simultaneously derive the 
kinematical center, systemic velocity, position angle of major axis (P.A.) and inclination 
of galactic disk [i), in addition to the rotational velocity (Kot)- Initial gue sses for the 



fitting are the o bserving center for kinemati cal center, values from BIMA SONG. ISofue et al. 
(119991 ) ■ or RC3 (Ide Vaucouleurs et al.lll99ll ) for P.A., i, and systemic velocity. For the BIMA 



data, however, inclination angles were derived to be about 90 degree for almost all galaxies, 
although the sample was selected to be nearly face-on. This discrepancy is presumably due 
to the dumpiness of CO dist ributions, and we u se values of i from the literatures (BIMA 
SONG. ISofue et al.l (Il999l ). or lDaigle et al.l (120061 )). Derived parameters from this procedure 
are listed in Table [3], with adopted values for P.A. and i. 

Fixing all the parameters except for V^ot; the rotation curve with statistical errors can be 
calculated from the velocity field. We could not derive meaningful kinematical parameters 
for 3 galaxies (NGC 3184, NGC 4535, and NGC 4736) out of 13 galaxies, however. For 
NGC 3184 and NGC 4736, the kinematical center could not be determined. For NGC 4535, 
the kinematical center and the systemic velocity could be determined, but the rotational 
velocities were not properly derived at some radii. This failure could be attributed to that 
CO data for NGC 4736 are rather noisy and those for the other 2 galaxies do not contain 
single-dish data, so that the fiux from extended components are resolved out. For the 
remaining 10 galaxies we derive rotational velocities, but some of them are limited in the 
range of radii or have larger errors in the outer regions. 



Given that RCs from CO are not sufficient for all the galaxies, we searched the literature 
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for RCs with a wider range of radii. ISofue et al.l (119991 ) compiled data from optical (Ha and 
[Nil] line) and radio (CO a nd Hi line) observat ions to present central-to-outer rotation curves 
for nearby spiral galaxies. iDaigle et al.l (120061 ) observed a part of the SINGS galaxies in the 
Ha line via Fabry-Perot spectrometry and obtained kinematical information and rotation 
curves in the same way as in this study. 

In Figure [3l available rotational velocities multiplied by the sine of the assumed incli- 
nation an gle (Knt sini ) are p lotted against radius. In the central regions, rotational veloc- 
it ies from Sofue et al.l ( 19991) are lar g er tha n those from the others, while in outer regions 



Sofue et al. 



Jl999h and IPaigle et all (l2006h give almost consistent results. This difference 
is principally d ue to that ISofue et al.l (119991 ) traces terminal velocities in a position-velocity 
diagram, while iDaigle et al.l (120061 ) and this study trace intensity-weighted velocities, and 
that these two velocities are known to be oft en substantia l ly dif ferent in the inner regions. 
RCs from CO data generally follow those from lDaigle et al.l (120061 ) but drop around the outer 
edge of the CO map. 

From these available RCs, we selected which to use considering errors and radial cover- 
age, in the following analysis. In principle, we try to avoid the use of CO data, since they are 
not reliable in the outer regions. Although not all the data are consistent each other, most 
of them are in good agreement where we measure the offsets. For some galaxies, however, a 
poor RC makes the application of offset method difficult, indicating the uncertainty of the 
results could be underestimated. We could not obtain a satisfactory RC for NGC 4535, so 
we decide to exclude this galaxy from the following analysis. Available RC data and the 
choice of RC are summarized in Table [2l 



3.3. Phase Diagram and Offsets 

The image analysis was started by checking the coordinates of CO and Ha data. If the 
coordinates of Ha images are incorrect or not available, foreground stars seen in the i?-band 
image with the same field of view as Ha were used to derive the coordinates. The Aladin 
Sky Atlas was used to obtain the coordinates of foreground stars and the task ccmap of IRAF 
was used to calculate and register the image coordinates. 

Projected images of galaxies in the (RA, Dec) coordinate were transformed into depro- 
jected images in the polar coordinate (radius r, azimuth 0), by the use of the task pgeom 
of AIPS with the dynamical center derived above and adopted values of P.A. and i listed 
in Table [31 For NGC 3184 and NGC 4736, whose dynamical centers could not be derived, 
the center of the CO map, which is almost the same as the NED position, was adopted. In 
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Figure HI we show CO contours on an Ha image in the (0, logr) coordinates for 12 galaxies 
whose RC are available. Images in this coordinate, called the phase diagram, are used to 
recognize spiral arms, since spiral arms are often expressed as Inr = tan(zpitch)0, where ipuch 
is the pitch angle of arms, and appear as a line in this diagram. 

To measure the offset, we first derived the average flux density in r ± |Ar for each 
azimuthal angle. The step in radius Ar was set to a value ranging from one-third to half of 
the CO beamsize. Then, peaks of CO and Ha flux belonging to a spiral arm were searched, 
and the offset was measured as the azimuthal separation of CO and Ha peaks. 



3.4. The n-e Plot and Fitting 

The measured offsets (6) were plotted against Q in Figure O The error bar in 6 cor- 
responds to the spatial resolution of CO data, which is the largest factor of uncertainty in 
6. The uncertain ty in Q is calculated from errors of the RCs. Since errors of RCs from 



Sofue et al.l (119991 ) were not available, we assumed an error of ±20 km s~^. The least 
method was used to fit the Q — 6 plot and to derive Qp and tsF with errors. The inclusion 
of the uncertainty in Q and the use of the least method for the fitting are major revisions 
from what we have presented in the previous paper. Results from the fitting procedure are 
presented in §11 



4. Result 



As described in the previous section, we show a plot of 6 against Q for 12 sample galaxies 
in Figure [H For several galaxies with a substantial number (about 20 or more) of offsets, 
offset distributions for each arm are shown by different symbols. We found that the Q — 6 
distribution varies with galaxies, and that for some galaxies, 6 can even show a negative 
dependence on Q, which is not explained by our picture or Figure [H From these plots and 
the CO and Ha images, we categorized the sample galaxies into 3 types as 



C: galaxies with clear offsets, which show clear spiral structures in CO and Ha, and the fit 
to their Q — 9 plot results in tgp > and AQp/Qp < 1, 

N: galaxies with no offsets, which show clear spiral structures, but their offsets in arm 
regions are almost zero, 

A: galaxies with ambiguous offsets, which cannot be categorized into the above two, mostly 
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because their Q — 6 plot show a negative correlation or large dispersion (i.e., tsF < 
or AQp/Qp > 1). 

For our sample of 12 galaxies, we have 5 in C, 2 in N, and 5 in A. For the 5 galaxies 
in category C, we derived meaningful values of flp and tsF by the fitting. Results of the 
fitting and the offset category are listed in Table HI 

Results and discussions for each galaxy are presented in §4.H and characteristics and 
differences in physical and observational parameters of the 3 categories are discussed in §4.21 
Discussions on the derived Qp and tsF are presented in §5.11 and 15. 2[ respectively. 



4.1. Description for individual galaxies 

4.I.I. Galaxies with Clear Offsets 

Here, we describe our results for "C" galaxies. The Q — 6 plot with the fitted line is 
shown in Figure [6l In Figure [TJ we plot Q and Q± k/2 as a function of radius to locate the 
kinetic resonances. We could derive the lower and upper limit of Rqr for 4 out of 5 "C" 
galaxies, while only the lower limit was derived for the remaining 1 galaxy, NGC 4303, whose 
uncertainty in Qp is as large as Qp itself. On the other hand, the insufficient quality of RC 
for some galaxies or at some radii hampers reliable derivation of n, which makes it difficult 
to derive the Lindblad resonances. We could thus locate the ILR for NGC 5194 only. 

Comparison with previous studies is also presented in this section. Values of Qp and 
Rqr from this and other studies are listed in Table O 



NGC 0628 (M 74) We were able to measure 33 offse ts between CO and Ha at r = 0.5-3.7 



kpc (0'.2 — 1'.7), while Q was calculated from the RC of iDaigle et al.l (120061 ) with a correction 
for the adopted inclination angle of 24 degree. We should note that even though the Q — 6 
plot shows a large dispersion, the fitted line has a positive gradient since the outlying data 
with large errorbars have lower weights in the x^-fitting. Removing such data points with 
AQ/Q > 1 did not substantially change the fitted result. As offset distributions for the 2 
arms appear to be consistent, we derive only one value for tsF and Qp. 

From the derived value of Qp and the RC, the corotation radius is calc ulated to be RcTt — 
2—4 k pc, corresponding to 0'.9 — 1'.9. With RcTt — 5 kpc or 1^8 l ocated by lCepa fc Beckman 
( Il990l ) from the arm/interarm ratio of SFE, iFathi et al.l (l2007al ) derived flp = km s~^ 



kpc ^ with the same RC which we used. We should mention that the discrepancy in Qp is 
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due to the difference in the assumed distance and inchnation angl e, and that the location 
of the corotation resonance is consistent within errors. However, lElme green et al.l (Il992l ) 
examined the 5-band morphology of this galaxy and derived Rqr ~ 2'. 4, which is larger 
than our result. 



NGC 4254 (M 99) For this 3-armed spiral galaxy, we were able to measure 18 offsets 
between CO and Ha for the arml, extending from east to south, and arm2, extending from 
west to north, at r = 2 — 6 kpc (0'.4 — 1'.3). The difference in the Q — 6 distribution for the 2 
arms is not clear, since the number of offsets is small for arm2 (see Figure E]). We could not 
measure any offsets for arm3, originating from the same position as the arml but with much 
smaller pitch angle extending to north, since the Hll regions were not well aligned where the 
CO emission delineates the arm. 

This galaxy has an asymmetric structure, as one spiral arm (arml) is much prominent 
than the other two, which could be explained by a ram pressure eff ect from intra-cluster 
matter in the Virgo cluster (IHidaka fc Sofud l2002l : ISofue et al.ll2003bl ). In addition, a faint 
but large Hi structur e, VIRGOHI 21, which might be interacting with this galaxy, has 
recently been found (IHaynes et al.l 120071 : iMinchin et al.l 120071 ) and is also thought to be 
responsible for the asymmetry. This global asymmetry could explain why only a few offsets 
could be measured in arm2, and none in arm3. 



R 



From the derived value of Qp and the RC, the corotation radius was calc ulated to be 
~ 4.5 — 6.0 kpc, corresponding to 1'.0-1'.3. Our result is consistent with iKranz et al. 



( l200ll ). who used results from hydrodynamic simulations to locate the -Rcr- 



In our previous work, f2p and tsF were derived to be 28^1^ km s^^ kpc^^ and 4.8 ± 1.2 
Myr, respectively, while this work gives 10 ± 3 km s^^ kpc~^ and 12.4 ± 1.3 Myr. The 
difference between the previous and current work is the fitting scheme, CO data (both in 
map and RC), and the adopted inclination angle. We have included the uncertainty in Q in 
the fitting procedure, and furthermore, the number of data and the radial range of RC have 
been increased by our new CO observations with 3 pointings with NMA. These two aspects 
should have improved the reliability of results. The assumed inclination angle is changed 
from 34° to 52.4°. However, since sin(34°)/ sin(52.4°) ~ 0.7, this is not enough to explain 
the discrepancy in Qp and tsp (See §5.31) . This implies that the uncertainties from the fitting 
could be underestimated. 



NGC 4303 (M 61) We were able to measure only 8 offsets between CO and Ha at 
r = 2.5 — 4 kpc (0'.5 — 0'.9), so that the uncertainty in flp is too large to derive Rcr for this 
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galaxy. We only give the lower limit of Rcr ^ 3 kpc or 0'.6. 



As NGC 4303 has a double bar and an active nuclei, studies o f this galaxy ha v e con- 
centrated on the dynamics and morphology in the central region. IColina &: Wadal (120001 ) 
performed multiphase hydrodynamical simulations and derived Qp ~ 0.5 Myr^^ ~ 500 km 
s~^ kpc~^ for the nuclear bar. They loc ate the CR at the end o f the nuclear bar, which is 
about 2" in radius. On the other hand, iRautiainen et al.l (120051 ) derived Rcr = l'-5 ± O'.l 
for the outer bar with the size of 1'.5. Although their Rcr gives flp close to our mean value 
of 24 kr n s~^ kpc~^, it i s still unclear whether the bar and spiral arms have the same pattern 



speed. iKoda fc Sofud (120061 ) calculated orbits in the bar fr om CO data, taken as a part 
of the Virgo High- Resolution CO Survey (ISofue et al.ll2003al ). and estimated ages of young 
clusters apparent in Ha to be < 10 Myr, which is consistent to our result of tsF — 10 Myr. 



NGC 5194 (M 51) We were able to measure 41 offsets between CO and Ha at r = 2.5 — 6 
kpc (0'.9 — 2'. 2), and found that their dependence on Q is different for the 2 spiral arms. 
The offsets of arm2, directly connected to the companion galaxy NGC 5195, showed negative 
dependence on Q. We have therefore derived flp and tsF by the x^-fitting to the offsets of 
arml only. The difference in the Q — 6 distribution for the 2 arms, may be attributed to the 
tidal interaction with the companion. This is because the interaction can cause asymmetry 
in the disk, and the outer part of the spiral arm connected to the companion could be a 
material arm. 

From the derived f2p and the RC, we calculate Rcr = 8.0^°!o ^P"^' '^'■^-ol^ the 
ILR to be at r ~ 2 — 4 kpc, which is slightly outside of the inner edge of spiral arms. The 
derived values of -Rc r in Table E] span a ra nge more than a factor of 2, indicating the difficulty 
of its determination. IZimmer et al.l (1200411 derived flp to be abou t 40 km s~^ kpc~^, from the 
TW method using the 45m CO data jNakai et al.lll994h. w hile iMeidt et all J2008h applied 
the radial TW method to other CO data (IShetty et al.ll2007l ) and derived two pattern speeds 
of about 90 and 50 km s~^ kpc^^ for r < 2 kpc and 2 < r < 4 kpc, respectively, with a 
possible third patter n of Qp ~ 20 km s~^ kpc~^ for 4 < r < 5 kpc. As the results from 



Zimnaer et al.l (120041 ) and this work are both between 50 and 20, de rived bv iMeidt et al. 



(120081 ) . these three works are in good agreement. On the other hand, iKnapen et al 



Jl992h 



calculated the arm/interarm ratio of SFE al ong the two arms , and showed that positions of 
ILR, CR (at r ~ 2' or 6 kpc), and OLR from lElmegreen et al.l (119891 ) were in good agreement 
with the minimum of this ratio. They also found that the pattern of this ratio along the 
position in arms is quite similar for the two arms, and claimed that some global triggering 
mechanisms for star formation should be at work. This pattern they found is not consistent 
with our result showing different offset distributions for the two arms, as well as the location 
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of the CR. 



Bastian et al.l (120051 ) estimated ages of Hll regions in this galaxy as < 10 Myr based on 



population synthesis models, and the derived value of tsF is consistent to their ages. Recently 
published extinction-corrected distribution of star forming regions ()Kennicutt et al.l 120071 ) 
and the latest CO observations with the 45m and CARMA (IKoda et al.ll2009l ) of this galaxy 
will be able to improve the current results in terms of reliability of location of Hll regions, 
the spatial resolution of CO map, and the quality of the rotation curve. 



NGC 5457 (M 101) We measured 10 offsets at r = 1.7 - 3.1 kpc (0'.8 - 1'.5). Although 
the data range in Q is rather narrow, we were able to derive tsF and Qp. We should note 
that the region where we measured offsets is only the central part of the whole disk, whose 

R25 is ii'.g. 



Cedres fc Cepal ((2002) derived the Ha equivalent width from observations with narrow 



band filters, and estimated a mean age of Hll regions at r < 5' to be about 1.6-4 Myr, 
supporting our result of tsF = 4.0 ±1.3 Myr. 

The fitted Qp of 72 ± 3 7 km s "^ kpc^^ locates the corotation at Rqr = 2.7lo;g ^P^^; 
1.31q 4 arcmin. IWaller et al.l (119971 ) derived f2p = 19 ± 5 km s~^ kpc"^ using offsets between 



CO and FUV at r ^ 1^5, assuming the timescale from CO to FUV to be 3 Myr, locating 
RcR = 5'.5±1'.5. lElmegreen et al.l (119921 ) also located the corotation at 5'. 5 from the optical 
appearance. This discrepancy implies a radial variance of the pattern speed or the inclination 
angle, which could be caused by interactions with companions, NGC 5474 and NGC 5477. 
If the inner disk where we measured offsets, is more inclined than the adopted inclination 
angle of 18°, the resultant Qp would be smaller or closer to other results mentioned above 
than the value presented here. In terms of the inclination, the adopted value for this galaxy 
is one of the smallest in our sample, implying a difficulty in deriving RC, which could be also 
inferred by the poor consistency in available RCs (see Figure E]). The uncertainties derived 
by the fitting, thus, could be underestimated. 



4.1.2. Galaxies with No Offsets 

Here, we describe our findings for "N" galaxies in comparison with previous researches, 
and discuss why they do not show offsets. Values for Qp from the literature are listed in 
Table El 
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NGC 4321 (M 100) We measured 34 offsets at r = 1.6 - 6.6 Icpc (0'.3 - 1'.4), and found 
that tfie offsets in the arm region (r ~ 5 — 6.5 kpc) are almost zero, while offsets close to the 



bar re gion (r < 4.5 kpc) are about 10 degrees, which had already been found by lSheth et al. 



( 120021 ). There are several possible reasons for no offsets found in the arms: (1) material arms, 
(2) corotation at where offsets were measured, and (3) elliptical orbits nearly parallel to the 
spiral arms. Since there are offsets close to the bar region, where orbits should have higher 
ellipticity than in the arm region, (3) seems to be less plausible at least for this galaxy. 



+4 1,^ ^-1 



Rand &: Wallin] (12004J ) derived the pattern speed of this galaxy to be 2 8_c, km 



kpc ^ with the TW method applied to the CO data. With the rotation curve from lSofue et al. 



( 119991 ) ■ this value locates the CR at about 10 kpc (~ 2') in radius. Previous studies listed 
in Table M also suggest Rcr ^ 2'. These measurements are well beyond where we measured 
offsets, so cannot explain the lack of offsets. 

The northern arm is weaker than the southern arm both in CO and Ha. This asymmetry 
is thought to be due to the existence of the companion galaxy NGC 4 322, the ram p r essure 



from intragalactic matter in the Virgo cluster, or the central bar. iKnapen et al.l (119961 ) 
investigated the arm/interarm ratio of the SFE, and found no symmetric patterns though 
this ratio was larger than unity for most part of the two arms. They concluded that this 
was because their observed region did not include any resonances, although star formation 
on arms was actually enhanced and it might be triggered by spiral density waves. Thus, 
we deduce that the spiral density wave might be rather unstable or localized at least in the 
region we see, and that (1) would be the most plausible among the three. If this grand-design 
galaxy is confirmed to be without density waves, it will give a new picture of spiral galaxy 
formation. 



NGC 5248 We were able to measure 13 offsets at r = 4.2 — 9 kpc (0'.6 — 1'.4), which were 
found to be almost zero, while 2 spiral arms were clearly traced. Three possible reasons for 
no observed offsets in NGC 4321 are also applicable to this galaxy. 



In addition to a small bar with size of 0'.4, iJogee et al.l (120021 ) found a large-scale stellar 
bar with a semimajor axis of 1'.6, indicating that the spiral structure seen in CO and Ha could 
be inside the bar. They also estimated -Rcr ^ 1^9 from the as s umed ratio of the corotation 
radius to the bar semimajor axis (1.2), while lElmegreen et aLl (Il992[) derived a comparable 



Yuan fc YanJj2006h developed a nonlinear 



value. Following the concept of the larger bar, 
model for the spiral density waves excited by a bar potential, and succeeded in explaining 
the observational features of this galaxy. Their orbits show strong inward streaming motions 
along the spiral arm at r < 0'.8, or 5.8 kpc using our distance. As we found that offsets of up 
to 1'.4 in radius were almost zero, this non-circular orbit can partially explain this feature. 
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4.1.3. Galaxies with Ambiguous Offsets 

In this section, we describe our findings for "A" galaxies in comparison witli previous 
studies, and discuss possible reasons for their ambiguity. These reasons are not the same for 
all the galaxies, and we expect that if data with higher quality become available, some of 
them will be recategorized as C. 



NGC 3184 We measured 19 offsets at r = 1.6 — 3.6 kpc (0'.6 — 1'.4), and found a negative 
correlation between 9 and with large dispersion. 

Since CO data from BIMA SONG do not include data from single-dish observations, 
they are subject to the missing flux problem. In order to quantify how much flux is missed 
i n the BIMA obse rvations, we compared its spectra to those from the 45m observations 



( iKuno et al.l 120071 ) in Figure [H While about 90% of the total flux was detected in the 



central 21" region, about 70% was missed if the flux was averaged at r < 45 

We do not, however, attribute the ambiguity in the offset to the missing flux, since spiral 
arms are thought to be comprised of smaller molecular clouds, which should be selectively 
detected in the interferometric observations. Since the spatial resolution of the CO data 
(~ 250 pc) is small enough to resolve typical molecular arms and offsets, the field of view 
and/or sensitivity should be insufficient. As spiral arms can be traced in the Ha to larger 
radii, we expect that the ~ 9 distribution will become clearer with smaller dispersion, if 
CO data with a larger field of view and higher sensitivity become available. 



NGC 3938 We measured 9 offsets at r = 3 — 4.6 kpc (0'.6 — 0'.9), and found a negative 
correlation between 9 and fi, with a large dispersion. This large dispersion is partly because 
of the large errors and unreliability of the CO rotation curve, w hich was used to calculate 
VL. The reason we did not use the RC from iDaigle et al.l (120061 ) is that it is systematically 
different from ours, suggesting that the inclination angle mentioned in their paper is not 
correct. In addition, the CO image and resolution around 500 pc indicates that molecular 
arms are not sufficiently resolved and detected, while spiral arms are clear in the Ha image. 
Thus we need further CO observations to discuss the spiral structure of this galaxy based 
on the VL — 9 distribution. 



NGC 4736 (M 94) Though we were able to measure 8 offsets at r = 1-1.5 kpc (0'.7-l'.0) 
which are clear in the CO and Ha images, they show a negative dependence on VL. We derived 



the ro tational velocity from CO data, but it was not consiste nt with the RC fro m lSofue et al. 



(Il999l ) where we measured offsets. Given that the RC from lSofue et al.l (119991 ) is flat, which 
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is typical of spiral galaxies, we use this to calculate f2. However, even if we used the RC 
from CO data, the dependence of offsets on Q would still be negative. 

NGC 4736 is the only galaxy in our sample to be classified as flocculent (AC=3). A star 
forming ring or short spiral arms are apparent at r ~ 1 — 1.5 kpc in the Ha image, while 
molecular arms are longer and apparent at r ~ 1 — 2 kpc. Th i s diff erence in morphology 
between CO and Ha, which has been also noted by lWong &: Blita (l2000l ). could be responsible 
for the negative correlation between 6 and Q. Therefore, we deduce that there are other 
triggers of star formation than the spiral density wave in this galaxy. 



One p ossible trigger is a n exp anding poststarburst ring, proposed by Ivan der Kruit 



(jl974j ) and ISanders &: Banial (119761 ) . They explained that the ring resulted from a central 
starburst which occu rred about 10 Myr ago, from spectroscopic observations and hydrody- 
namical calculations. Maoz et al. found two compact UV sources separated by 2". 5 in 

the nucleus with diffuse emission centered on the one source, presumably corresponding to 
the optical nucleus. They interpreted their results to be the final stag e of a merger, which is 
consistent with the poststarburst scenario. iMuiioz-Tuiion et al.l (120041 ) adopted another pos- 
sibility of this ring to be located at the OLR of the central bar as well as the ILR of the outer 
disk. The pattern speed for each structure should be 85 and 27 km s~^ kpc~^, respectively, 
and the latter could explain the location of outer faint ring at r ~ 5'. They also noted that 
the FUV ring is located slightly outward from the Ha ring, implying an inward-propagating 
star formation, which is totally opposite with the poststarburst scenario. 



Rand fc WallinI (120041 1 applied the TW method to BIMA CO data, and derived 



152 ± 28 km s~^ kpc~^, correc t ed fo r our distance. This value is much larger than that 
derived by iMufioz-Tufion et al.l (120041 ). and locates the corotation at r = 0'.8 — 1'.2, which 
is almost same as where we measure the offsets. 



NGC 6181 We were able to measure 27 offsets for 2 spiral arms at r = 2—5 kpc (0'.2— 0'.5). 
While the offsets of arml show a negative dependence on Q, those of arm2 show a positive 
dependence (Figure [5]). We thus fit the Q — 9 plot for arm2 only, but the uncertainties of 
the resultant values for tsF and Qp are too large to give a quantitative conclusion. As the 
large uncertainty is mainly due to large errors in the RC we used, further observations will 
be able to give more reliable results at least for arm2. 



NGC 6946 At the region where CO emission was detected, the spiral structure is not 
clear, so we were able to measure only 7 offsets at r = 1.1 — 2.4 kpc (0'.7 — 1'.5). This small 
number of data resulted in the large uncertainty of flp, so we categorized this galaxy as A. 
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The positive slope in the — 6 plot, however, gives the fitted value of tsF = 1-1 ± 0.8 Myr. 
This is slightly smaller than that for C galaxies, but still has the same order of magnitude. 



Zimmer et all (120041 ) applied the TW method to the IRAM 30m CO data flWalsh et al. 
2OO2I ). and derived r ^p = 39 ± 8 km s^^ kpc~^, consistent with our mean value of 38 km 



kpc \ iFathi et al.l (l2007bl ) also applied the TW method using Ha data from iDaigle et al. 



(l2006l ). whose spatial resolution was higher than that of CO data, and obtained two pattern 
speeds: Qp = 4:7^2 km s~^ kpc~^ for the inner structure (r < 1') and Qp = 22^^ km s~^ 
kpc~^ for the outer structure (r > 1' ). They also claimed that spatially smoothed Ha data 
gave a comparable value to flp from IZimmer et al.l (120041 ). As we measured offsets around 
the border of the inner and outer structure, our Qp, which is in between the two values, is 
not inconsistent to their results despite the large uncertainty. 

Since the spir al structure is mor e conspicuous and stronger CO emission is detected in 
the outer regions (IWalsh et al.l 120021 ). further CO observations pointed to such regions will 
give understanding of the property of the spiral arms and star formation in this galaxy. We 
anticipate that this galaxy will be recategorized as C if we can measure a larger number of 
offsets and thus derive Qp with smaller uncertainty. 



4.2. Property of Galaxies in Each Category 

In order to understand relationships between spiral structures and other physical param- 
eters, we have examined several properties of galaxies according to their category. Similar 
analyses are performed for all "C" galaxies, and their results are shown in §5.2[ 



4.2.1. Molecular Gas Property 

As we found that the missing flux in CO data of NGC 3184 data could be as large as 
about 70%, we excluded this galaxy from the following quantitative analysis of molecular 
gas. NGC 6181, whose CO data are also from interferometric observations only, is included 
in th e analysis, since we find almost no flux is missed in the data by comparison with spectra 



from lYoung et al.l (119951 ) . 



First, we reassess the quality of CO data, which can cause biases to the categorization. 
We list the beamsizes (bmaj) and the noise levels (rms) of the CO images together with the 
offset category in Table [TJ and show a plot of rms against bmaj in Figure [Hi From this flgure, 
we flnd no systematic differences or biases to the category. 
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Next, we examine global properties for each category. We calculate the mean H2 surface 
density S using a conversion factor of Xqq = 3. Ox 10^*^ cm~^ (K km s~^)^^ at the radial region 
where we measured offsets. In addition, we derive S only for the arm regions, centered on 
CO peaks, which were used to derive offsets, with a width of 500 pc. From Figure [101 where 
disk-averaged and arm-only E are plotted for each category, we find no clear dependence of 
both values on the offset category. This indicates that larger (> 500 pc) scale properties of 
gas are not closely relevant to the spiral structure. 



4.2.2. Morphological Property 

In order to quantify the morphological properties of the underlying gravitational poten- 
tial of the spiral disk, we used the arm/interarm ratio in K-hand images as an indicator of 
the strength of the galactic shock or spiral arms, since the luminosity in K is less sensitive to 
extinction and thought to be proportional to the total mass of low-mass stars, which domi- 
nate total stell ar mass, assura i ng a constant mass-to-luminosity ratio. The K-hand images 



are taken from iKnapen et al.l (120031 ) and 2MASS, and shown in Figure [TTl We defined the 
arm region as a region with 500 pc width centered on the Ha peak which was measured 
in order to derive offsets, while the interarm region was set to the same width centered at 
the midpoint of the two arm regions, or 90 degrees away from an arm region if only one 
arm region was defined at that radius. The arm/interarm K flux ratio was calculated at 
each radius where Ha peaks (i.e., offsets) were measured and averaged over the disk. In 
Figure [121 the averaged ratio is plotted for each offset category. Although there is no strong 
correlation, "A" galaxies tend to have a smaller ratio. 

As an alternative way of estimating the arm strength, we performed an one dimensional 
discrete Fourier transformation, 

360-1 

A{r,m) = F(r, 0) exp(— 27rim0/36O), (7) 
4>=o 

where F{r, (p) is the i^'-band flux at the polar coordinate (r, 0), in order to know which mode 
m is dominant within a galaxy and obtain the typical strength of asymmetric components. 
We calculate the relative amplitude \A{r,m)\/\A{r,m = 0)| for m = 1 — 4, and the mean of 
these components is shown against the offset category in Figure [T3l The m = 2 component is 
the strongest for most of the galaxies, regardless of category. We find no clear dependence of 
the amplitude on the category, similar to that in the K flux ratio. The smallest amplitudes, 
however, belong to "A" galaxies. 



This tendency between the arm strengths and the offset category implies that stronger 
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stellar arms or a deeper potential generates clearer offsets. As the potential depth should be 
tightly correlated with the strength of spiral density waves, the ambiguity in offsets for "A" 
galaxies at least partially results from weaker density waves. 

In Table [HI we list the number of SA, SAB, and SB galaxies in each category. While our 
sample is not large enough for a robust statistical study, "N" galaxies seem to have stronger 
bar than others. Thus, bars could be responsible for no offsets as discussed in §4.1.21 



5. Discussion 

For 5 galaxies in the "C" category, we were able to derive their Qp and tsF by the 
X^-fitting. In this section, we discuss the dependency of the derived values on other physical 
parameters. As we have only 5 galaxies, the statements here are still tentative and require 
a larger sample to confirm them statistically. 



5.1. Pattern Speed 

In Figure [TH we plot ellipses on CO and Ha maps to show the location of the derived 
corotation, except for NGC 4303, whose uncertainty in Qp is too large to derive a reliable 
RcR- We find that -Rcr is nearly consistent with the edge of CO data, while Ha arms are 
more extended. This is a natural result from the fact that we did not find large negative 
offsets, but we cannot exclude the possibihty that the small CO field of view has caused a 
bias to this result. For a more robust determination, larger CO maps are needed. In Table 



ini the radius of the optical disk R25 from Nearby Galaxies Catalog (iTuUyl Il988l ) and the 
derived -Rcr are listed. From the ratio Rcr/R25, we find that the corotation radius is about 
half of R25, except for NGC 5457. We should note that this result is fully independent of 
the assumed distance, and less sensitive to the assumed inchnation angle. 



5.2. Star Formation Timescale 

As the derived values of tsF are about 5 — 30 Myr and the free-fall time (~ l/y/Gp) of 
a typical molecular cloud (p ~ 100 cm~^) is about 6 Myr, we deduce that the gravitational 
instability should be a dominant trigger for star formation in the spiral arms. Here, we ex- 
amine how the gravitational instability and other factors affect the star formation timescale. 
Similar analyses are performed for each offset category and their results are shown in §4.21 
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We should bare in mind that tsF is hnearly dependent on the assumed distance, and 
also dependent on the assumed inclination, where its error is substantial for galaxies with 
low inclination angles. Its dependency is discussed in §5.31 



5.2.1. Gravitational Instability 

In Figure [151 the mean H2 surface density S derived from CO data is plotted against 
tsF- Values in the left panel are derived by averaging the data azimuthally, while those in 
the right panel are derived only from arm regions. As we obtained tsp consistent with the 
gravitational collapse timescale, we expect larger density with smaller tsF- However, there is 
no clear correlation between tgp and S, suggesting that molecular properties at scale > 500 
pc (adopted width of arm regions) are not directly relevant to star formation activities. Since 
the typical size of a GMC is about 50 pc, which is much smaller than the spatial resolution 
of CO data we used, we need CO data with higher resolution as well as larger sample of 
spiral galaxies for further study. 



5.2.2. Metallicity 



We plot the heavy element abundance in the form of 12 + log(0/H) from lZaritsky et al. 



( I1994J ) against tgp in Figure [T6l Each galaxy has 2 metallicity values measured at r = 0.4i?25 
and r = O.SRs, where Rg is the scale length. We measured offsets around 0.8-Rs with the 
exception of NGC 5457, for which we used the region of r < 1'.5 while O.SRg = 1'.7. Since 
heavy elements can help the cooling and collapse of molecular clouds by their transitional 
lines, we expect smaller tgp with higher metallicity. Such a correlation is, however, not 
apparent here, implying that these effects by heavy elements are not critical to star formation 
processes. 



5.2.3. Spiral Strength 



In Fi gure [T7I and [T8l we plot the arm/interarm ratio and the Fourier amplitudes \ A{r, m)\/\A{r, m 
0)1 for m = 1 — 4 from i^'-band images against tsF- Neither plots show strong dependence 
on tsF- If we can confirm that the spiral arm strength has no correlation with tgp, the spiral 
shock or potential should have no direct effect on star formation and its role could be just to 
agglomerate small pre-existing molecular clouds t o form massive r nolecu lar clouds, in which 
star formation eventually occurs, as suggested by iKoda fc Sofud (120061 ). This is important 
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for understanding the relationship between galactic dynamics and star formation activities, 
and we need a larger sample of galaxies for further study. 



5.3. Validity of Offset Method 

Out of 12 nearby spiral galaxies, we have succeeded in deriving Qp and tsF for 5 galaxies, 
which are referred as "C" galaxies. The rest of the sample, comprised of 7 galaxies, are 
categorized as "N" or "A" according to their offset distributions. 

We have examined the surface mass density S of the molecular disk, and strength of 
spiral arms in §4.21 to find any differences in physical properties between the categories. 
While no clear difference is found in the molecular properties, there is an indication that 
galaxies which show ambiguous offsets tend to have a smaller stellar density contrast between 
arm and interarm regions. Thus, the i^-band morphology can be used as an additional 
selection criterion, when a larger sample is needed for further study. 

The validity of the assumptions we have used in the offset method (i.e., circular rotation, 
constant tsF and Qp) can be estimated from the Q — 6 plot. Several features found in this 
work, such as the lack of offsets and different distributions for individual arms, indicate 
that these assumptions or even the density wave theory are not applicable to all the spira. 1 



galaxies. A presence of a central bar (e.g., lAthanassoula 



1980 



Romero-Gomez et al.l 120071 ) 



and an interaction with companion galaxies (e.g.. lToomre fc Toomrdll972l ) can drive material 
or short-lived spiral structures or enhance pre-existing density waves. Possible explanations 
for individual galaxies are mentioned in §4.11 but a full consideration of these effects is 
beyond the scope of this paper. 

Here we consider the effects of the adopted parameters on the results. The parameter 
we must care most is the inclination angle i, as it changes both tgp and flp. Since Q is 
derived as Kbs/ sini/r, larger tgp and smaller Qp will be derived, if larger i is adopted. In 
addition, different i can change the shape of spiral arms, since r = robs/ sini for the minor 
axis direction. This change is difficult to quantify, because it does not change anything 
regarding the major axis direction. These effects are especially significant for the most face- 
on (small i) galaxies. For example, if i = 12° is adopted for NGC 0628, which is half of 
our value, the resultant Qp and tsF will be changed by a factor of 2 from our result, since 
sin(12°)/ sin(24°) ~ 0.51. While this can directly change the result on tsF, the location of 
the corotation is almost independent of i, since fl is also inversely dependent on sini as well 
as Qp. 



Another parameter we must keep in mind is the distance D. As r is proportional to D, 
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Q is inversely proportional to D. Though D affects and tsF in the same way as sini, its 
uncertainty is about 20%, so that the effect is not as large as sinz for small i. We should 
note that the location of corotation is fully independent on D. 

In addition, the poor determination of the RC for some galaxies and the discrepancy 
between this work and previous work for NGC 4254, both indicate that the uncertainty of 
the fitted Qp and tsF could be underestimated. The latter and several other works also 
suggest a radial variation of Qp. In order to solve these problems, an extensive CO survey 
of nearby spiral galaxies with both interferometers and single dishes is critical. 



6. Conclusion 



We have presented a revised method to determine both the pattern speed (flp) and sta r 



formation timescale (tsp) of spiral galaxies simultaneously, proposed by lEgusa et al.l ( 12004| ). 
This method utilizes offsets between molecular and young-stellar arms, and we refer to it as 
the "Offset Method". 

We applied the offset method to CO and Ha data of 13 nearby galaxies, and derived 
Qp and tsF for 5 galaxies. Since their offsets are clear, we categorize them as "C" galaxies. 
From the results of these 5 galaxies, we find the followings. 



• The corotation radius calculated by the derived fip is near the edge of CO data, and 
is about half of the optical radius for 3 galaxies. 

• The derived tsF is roughly consistent with the typical free-fall timescale of molecular 
clouds, which indicates that gravitational instability is a dominant mechanism trigger- 
ing star formations in spiral arms. 

• The surface density of molecular gas calculated from CO data, heavy elements abun- 
dance, and spiral arm strengths evaluated from i^-band images do not show clear 
dependence on tsp- 



We also find that 2 out of the remaining 8 galaxies show no offsets between CO and 
Ha and categorized them as "N" galaxies. As their arms are clearly traced, the spiral density 
wave is usually thought to be at work in these galaxies. There are several possible reasons 
to explain the lack of offset in these galaxies: (1) material arms, (2) corotation where arms 
are seen, and (3) elliptical orbits nearly parallel to spiral arms. Although the current data 
are not sufficient to confirm which case is at work, we have gained an insight that a central 
bar could account for this feature, since these 2 galaxies are both barred. 
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With 1 galaxy excluded from our analysis due to its poor rotation curve, the remaining 
5 galaxies have ambiguous offsets, whose dependence on the rotational frequency is not clear 
or different from expectation. We categorize them as "A" galaxies. The major reasons for 
this ambiguity are (1) the density wave is weaker, and/or (2) observational resolution and 
sensitivity are not sufficient to detect arms and their offsets clearly. The former is supported 
by our finding that the arm strengths of "A" galaxies are slightly weaker than those of "C" 
galaxies. 

In addition to the results for individual galaxies, properties of molecular gas and mor- 
phology were examined according to the category. We list the results in the following. 

• The disk- averaged and arm-only value of molecular surface density do not appear to 
correlate with the category, indicating that larger (> 500 pc) scale properties of gas 
do not vary according to a spiral structure. 

• Analyses of the K-hand images show that the amplitude of stellar spiral arms is slightly 
smaller for "A" galaxies. 

Prom the results for individual galaxies and offset categories as mentioned above, we 
summarize our findings as 

1. Star formation in spiral arms is dominantly triggered by gravitational instability of 
molecular clouds. 

2. Stellar spiral arms could affect the appearance of offsets between molecular gas and 
young stars. 

We are grateful to Dr. Rebecca Koopmann, Dr. Crystal Martin, Dr. Johan Knapen, 
and Mr. Olivier Daigle for kindly providing their data and advices on how to use them. 
This work has made use of the NASA/IPAC Extragalactic Database (NED) and 2MASS 
Image Services via the NASA/IPAC Infrared Science Archive, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with the National 
Aeronautics and Space Administration; the VizieR catalogue access tool and the Aladin Sky 
Atlas, CDS, Strasbourg, France. F. E. and S. K. are thankful to the financial support by 
the Japan Society for the Promotion of Science during this research. A part of this study 
was also financially supported by the MEXT Grant-in- Aid for Scientific Research on Priority 
Areas No. 15071202. 
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Table 1. General Property of Sample Galaxies 



NGC 


R.A. (J2000) 


Dec. (J2000) 


Velocity 


Distance 


Morphology 


Tt 


AC* 




(h m s) 


(d m s) 


(km s-i) 


(Mpc) 








0628 


01 36 41.77 


+ 15 47 00.5 


657^ 


7.3^ 


SA(s)c 


5 


9 


3184 


10 18 16.98 


+41 25 27.8 


592=* 


8.7^ 


SAB(rs)cd 


6 


9 


3938 


11 52 49.45 


+44 07 14.6 


809^^ 


17^ 


SA(s)c 


5 


9 


4254 


12 18 49.63 


+ 14 24 59.4 


2407^= 


16.1'= 


SA(s)c 


5 


9 


4303 


12 21 54.90 


+04 28 25.1 


1566'' 


16. 1^= 


SAB(rs)bc 


4 


9 


4321 


12 22 54.90 


+ 15 49 20.6 


157P 


16. 1^= 


SAB(s)bc 


4 


12 


4535 


12 34 20.31 


+08 11 51.9 


1961^ 


16. 1^= 


SAB(s)c 


5 


9 


4736 


12 50 53.06 


+41 07 13.7 


308"^ 


5. 1'' 


SA(r)ab 


2 


3 


5194 


13 29 52.71 


+47 11 42.6 


463^ 


9.6'^ 


SA(s)bc 


4 


12 


5248 


13 37 32.07 


+08 53 06.2 


1153^ 


22.7^ 


SB(rs)bc 


4 


12 


5457 


14 03 12.59 


+54 20 56.7 


24P 


7.2'^ 


SAB(rs)cd 


6 


9 


6181 


16 32 20.96 


+ 19 49 35.6 


2375^ 


36.7^ 


SA(rs)c 


5 


12 


6946 


20 34 52.34 


+60 09 14.2 


48^ 


5.5'^ 


SAB(rs)cd 


6 


9 



Note. — (RA, Dec) and morphology are from NED. f: T is a numeric index correspond- 
ing to the morphology, and i s called the Hubble T type. |: AC is the arm class from 



Elmegreen fc ElmegreenI (|l987f ) 



Reference s. — (a) iHelfer et all |2003l ). (b) ISofue et al.l (|l999l ). (c) 
Tullvl (Il988h 



Sofue et al. 



(2003a|), (d) 



Table 2. Property of Data for Sample Galaxies 



NGC 






CO 






Ha 




RC* 






Ullilll 


Y"yTn c 1 

1 ilib 


LCI. 




tnl 


ref 






(") 


(") 


(mjy bcam^^) 




(") 








0628 


7.2 


5.3 


51 


BIMA+12m 


0.43 


CTIO 


1 


6, CO 


3184 


5.9 


5.4 


50 


BIMA 


1.7 


JKT 


2 


6, CO 


3938 


5.9 


5.4 


59 


BIMA+12m 


0.3 


KP2 


1 


6, CO 


4254 


4.8 


3.6 


15 


NMA+45m 


2.1 


KP9 


3 


CO 


4303 


7.3 


5.5 


47 


BIMA+12m 


1.4 


JKT 


2 


7, CO 


4321 


7.2 


4.9 


50 


BIMA+12m 




S90 


4 


6, 7, CO 


4535 


7.3 


5.7 


63 


BIMA 


1.9 


JKT 


2 


CO 


4736 


6.9 


5.0 


64 


BIMA+12m 


1.4 


JKT 


2 


7, CO 


5194 


5.8 


5.1 


61 


BIMA+12m 




BS 


4 


6, 7, CO 


5248 


6.9 


5.8 


38 


BIMA+12m 


1.4 


JKT 


2 


CO 


5457 


5.7 


5.4 


45 


BIMA+12m 


2.4 


INT 


2 


7, CO 


6181 


3.4 


2.8 


11.5 


NMA 


4.7 


P60 


5 


CO 


6946 


6.0 


5.0 


61 


BIMA+12m 


1.4 


JKT 


2 


6, 7, CO 



Note. — f: The rms of CO data was calculated for a velocity width of 10 km s'-^. *: "CO" 
means a rotation curve from the velocity field of CO data, and data in bold font were used 
for the appHcation. (See N3.2l for detail.) 

Abbreviations for telescopes: BIMA (Berkeley-IUinois-Maryland Association), 12m (NRAO 
12m telescope), NMA (Nobeyama Millimeter Array), 45m (Nobeyama 45m Telescope), CTIO 
(Cerro Tololo Inter- American Observatory 1.5m), JKT (Im Jacobus Kapteyn Telescope), KP2 
(Kitt Peak 2.1m), KP9 (Kitt Peak 0.9m), S90 (Steward 90inch), BS (Kitt Peak Burrell- 
Schmidt), INT (2.5m Isaac Newton Telescope), P60 (Palomar 60inch) 



References. — (1) 


Kennicutt et al. 


(2003). (2) Knaoenetal. 


KooDmann. Kennev. & Young 


( 


2001). (4) 


Martin & Kennicutl 


(2001 


), 


et al.. private communication. f6) Daigle et al. ( 


20061). (7) Sofue et al. ( 


1999) 



( 20041) ■ (3) 
) Koopmann 
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Table 3. Parameters Derived from CO Velocity Field 



Dynamical Center Adopted Value 

NGC R.A. (J2000) Dec. (J2000) error Velocity P.A.t it P.A. i 

(hms) (dms) (km s"!) (°) (°) (°) (°) 



0628 


01 


36 


41 


589 


+ 15 


47 


01 


30 


3 


3" 


650.2+ 1.5 


11.8 + 1.1 




11.8 


24^ 


3184 




























176.4^= 


16. 7^= 


3938 


11 


52 


49 


502 


+44 


07 


20 


50 


9 


0" 


812.1 + 3.7 


202.2 + 5.1 




202.2 


24=^ 


4254 


12 


18 


49 


799 


+ 14 


25 


03 


90 


3 


2" 


2393.7+ 8.5 


70.9 + 4.1 


52.4+2.3 


70.9 


52.4 


4303 


12 


21 


54 


929 


+04 


28 


21 


60 


6 


1" 


1557.8 + 5.4 


304.8 + 9.1 


27.7+9.5 


318"^ 


27*^ 


4321 


12 


22 


55 


192 


+ 15 


49 


20 


10 


5 


7" 


1575.2+ 7.8 


161.7+9.5 


34.5 + 23.3 


146*^ 


27*^ 


4535 


12 


34 


20 


427 


+08 


11 


54 


10 


3 


2" 


1943.8 + 9.8 










4736 




























108'' 


35*^ 


5194 


13 


29 


52 


565 


+47 


11 


42 


20 


4 


2" 


469.2 + 2.4 






22^ 


20'' 


5248 


13 


37 


32 


140 


+08 


53 


03 


90 


7 


1" 


1165.0+10.9 


109.6 + 0.2 




109.6 


43^ 


5457 


14 


03 


12 


330 


+54 


20 


56 


99 


4 


3" 


258.1 + 2.3 






38^ 


18*^ 


6181 


16 


32 


20 


916 


+ 19 


49 


35 


23 





6" 


2395.7 + 4.0 


332.8 + 0.7 


66.8 + 3.9 


332.8 


66.8 


6946 


20 


34 


52 


703 


+60 


09 


13 


39 


4 


3" 


51.3 + 9.0 


244.8 + 15.7 






30^ 



Note. — f: P.A. is the position angle of the disk major axis, which is defined as 0° at north and increases 
counterclockwise, p. i is the inclination angle of the disk, which is defined as 0° when the disk is face-on, 
and 90° when edge-on. Data fields shown as " • • • " indicate that meaningful values are not derived by the 
analysis described in 



References. — (a) iHelfer et al.l (|2003l ). (b) ISofue et all (|l999l ). (c) iDaigle et al.l (|2006l ) 
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Table 4. Result of Fitting with Offset Category 



NGC 


tsF (Myr) 


rip (km s ^ kpc ^) 


^ of data 




Category'!' 


0628 


28.2 ±3.1 


16 ±3 


33 


91.78 


C 


3184 


-34.1 ± 10.9 


51 ±24 


19 


26.93 


A 


3938 


-6.4±6.5 


20 ±28 


9 


5.2 


A 


4254 


12.4± 1.3 


10 ±3 


18 


29.51 


C 


4303 


10.8 ±5.7 


24 ±29 


8 


7.09 


C 


4321 


2.3 ±0.8 


31 ±20 


34 


86.80 


N 


4736 


-17.2 ±4.8 


166 ±66 


8 


15.67 


A 


5194 


13.8 ±0.7 


40 ±4 


41 


811.54 


C 


(arml 


7.1 ±0.5 


31 ±5 


24 


157.32 


C) 


5248 


-0.1 ±2.6 


270 ± 10860 


13 


4.16 


N 


5457 


4.0 ± 1.3 


72 ±37 


10 


16.09 


C 


6181 


-85.7 ± 102.7 


38 ± 150 


27 


31.4 


A 


(arm2 


43.4 ±54.5 


33 ± 119 


14 


19.52 


A) 


6946 


1.1 ±0.8 


36 ±79 


7 


38.07 


A 



Note. — f : "C" galaxies are categorized if their tsF and ftp can be determined 
by the fl — 6 plot. "N" galaxies show no offsets between CO and Ha arms. "A" 
galaxies are those which do not fall into the above two categories. See ^for details. 
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Table 5. f2p Values for "C" Galaxies 



NGC Authority D P.A. i Hp _RcR Method Data 

(Mpc) {°) {°) (kms-lkpc-i) (arcmin) 

0628 Elmegreen et al. (19921 ■ ■ ■ ■ ■ ■ 21=^ ■ ■ • 2.4 morphology B 

Fathi et al. ('2007a') 9.7 25 6.5 31^3 l.S*" SFE Hi, Ha 

This work (2009) 7.3 11.8 24 16 ifc 3 0.9 -1.9 offset CO, Ha 

4254 Elmegreen et al. ( 1992) ~ ~ 27^^ ~ L5 morphology S 

Kranz et al. C2001') 20 67.5 41.2 20 1.3 ±0.2 simulation K,'Ra 

This work (2009) 16.1 70.9 52.4 10 ± 3 1.0-1.3 offset CO, Ho 

4303 Colina fc Wada (20001 16l 40 ~ 500 ¥' simulation 

^autiaine^_et_aL (2005) • • ■ ■ ■ ■ ■ ■ ■ ■ ■ • 1.5 simulation H 

This work (2009) 16.1 318 27 24 ± 29 > 0.6 offset CO, Ha 

5194 TuUy (1974) 4 ^10 20 90 2^2 morphology Ho 

Elmeg reen et al. (1992 ') ■ • ■ ■ ■ • 45^ ■ ■ • 2.2 morphology B 

Garci'a-Burillo et al. (1993') 9.6 170 20 27 2.7 simulation CO 

^oge^_et_aL (^1993) ■■■ -10 20 •■• 2.1 streaming motion Ha 

Kuno et al. (1995') 9.6 -10 20 14 5.0 morphology optical 

Oev et al. (2003) 9.6 ■ ■ ■ 5.0 isochrones Ha 

Zimmer et al. (2004*) 9.5 -10 20 38 ± 7 ■■• TW CO 

This work (2009) 9.6 22 20 31 ± 5 2.9lo 4 ^O, Ha 

5457 Elmegreen et al. ( 1992) ~- ~ 12^^ ~ 5l morphology B 

•Waller et al. (1997') 7.4 39 18 19 ± 5 5.5 ± 1.5 offset CO, FUV 

This work (2009) 7.2 38 18 72 ± 37 l-SlJi oS^et CO, Ha 



''Inclination is calculated from the ratio of major axis to minor axis, from lde Vaucouleurs et"al ] lfl97i 
''^CR is from lCepa &: BeckmanI l|l99Clh but Sip is derived from their new rotation curve. 



Table 6. Vt-p Values for "N" Galaxies 



NGC 


Authority 


D 


P.A. 


i 




RcR 


Method 


Data 






(Mpc) 


(°) 


n 


(km s^-"- kpc^-"-) 


(arcmin) 






4321 


Elmegreen et al. (1992') 






27=" 




2.0 


morphology 


B 




Sempere et al. (1995) 


20 


155 


32 


20 


1.8 


Canzian, simulation 


Hi 




Canzian & Allen (1997) 


17.1 


153 


28 


31 


1.6 ±0.2 


Canzian 


Ha 




"Wada et al. ( 1998) 


17.1 


155 


30 


65 (bar) 


0.2 (bar) 


simulation 


CO 




Oev et al. (2003) 


16.1 








2.6 


isochrones 


Ha 




Rand & Wallin (2004) 


16.1 


153 


27 


28+4 




TW 


CO 


5248 


Elmegreen et ah (1992) 










1.7 


morphology 


B 




Jogee et al. (2002) 


15 


135 






1.9 


bar length 


K 



''Inclination is calculated from the ratio of major axis to minor axis from lde Vaucouleurs et"ar 
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Table 7. Spatial Resolution and Sensitivity of CO data 



NGC 


bmajf (pcj 


rms* [Mq pc 


Category 


0628 


250 


52 


C 


3938 


490 


73 


A 


4254 


370 


34 


C 


4303 


570 


45 


C 


4321 


560 


55 


N 


4736 


170 


72 


A 


5194 


270 


80 


C 


5248 


760 


37 


N 


5457 


200 


57 


C 


6181 


600 


47 


A 


6946 


160 


80 


A 



^ Major axis of the synthesized beam from the CO 
data. 

■I- Sensitivity of CO channel maps (10 km s~^ per chan- 
nel) measured as the root mean square of emission-free 
regions. 

Table 8. Bar Classification and Offset Category 







N 


/V 


total 


SA 


3 





3 


6 


SAB 


2 


1 


2 


5 


SB 





1 





1 


total 


5 


2 


5 


12 



Table 9. R25 and Rcr 



NGC i?25 (arcmin) i?cR (arcmin) Rck/Ri5 



0628 


3.5 


~ 0.9-1.9 


0.3-0.5 


4254 


2.5 


~ 1.0-1.3 


0.4-0.5 


5194 


6.8 


2 9+°-2 


0.4-0.5 


5457 


11.9 


■'-•'5-0.4 


0.1 - 0.3 
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Fig. 1. — Basic idea of our method. If we observe a face-on spiral galaxy at t = (the left 
panel), the same galaxy will be observed as the right panel at tsF- The thick solid lines are 
molecular arms at time t of each panel. The thick dashed lines in the right panel (^sf) show 
the position of molecular arms in the left panel {t = 0). The offset distance between the 
massive stars and molecular arm is d, expressed in equation (jl]). 
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Fig. 2. — CO contours on an Ha image of sample galaxies. 




Figure [2]-continued 
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Fig. 3. — Available rotation cur ves from ISofue et al.l (119991 ). iDaigle et al.l (120061 ). and GAL 
applied to CO velocity field of Heifer et al.l ( 20031 ) or from our observations, labelled as 
"RC99", "SINGS", "BIMA", and "NMA", respectively. Labels in parentheses mean that 
their data are not available. 
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A^imutri [cogj A/'muin [acq] 



0, 0.4)1 O.0? D.Oi D.O' D.OS 0.08 0^7 K!0 BK) BOD lOCO i?M 'WO 




fl ISO .-ifiii 'i-'Lj iia &jo /^o 90 iao 2 /o j60 ^io Ma &io 720 



Fig. 4. — Phase diagram of each galaxy with CO contours on an Ho; image. The ordinate 
is the radius in logarithmic scale and the abscissa is the azimuthal angle and shown for two 
periods, or — 720 degree, to delineate spiral arms. 
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Figure Hl-continued 
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Fig. 5. — Plot of 6 against f2. Horizontal lines indicate 9 = 0. For several galaxies, two 
symbols (open and filled circles) are used to show difference between arms. 
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Fig. 6. — Q — 9 plot with the fitted fine for each galaxy in the C category. The gradient and 
horizontal-axis intercept of the fitted line correspond to ^sf and Qp, respectively. 
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Fig. 7 — Plot offl (solid line with '+'), + «;/2 (dashed line with 'x'), and fl-K/2 (dashed 
line with '*') as a function of radius. Derived Up is presented as a thick horizontal hne 
with thin solid lines indicating the uncertainty. The corotation and Lindblad resonances are 
derived to be where Qp = and flp — ± k/2, respectively. 
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NGC 3184 (center) 



NGC 3184 (r<45") 
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Fig. 8. — Spectra from BIMA data (dotted) and 45m data (solid). Flux is averaged for the 
central 21" region (left) and for the central 90" region (right). The missing flux derived by 
"1 - BIMA/45m" is larger when averaged for most of the disk than for the center. 
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Fig. 9. — Plot of sensitivity (rms) against beam major axis (bmaj) of the CO images listed 
in Table [7] for each category. 
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1 2 

Offset Category 




Fig. 10. — Mean H2 surface density against the offset category (1=C, 2=N, 3=A). Values in 
the left panel are derived by averaging azimuthally, while those in the right are derived by 
averaging only the arm regions. 



-43- 




W JB" J3" 




2t^ 11' rcr a*^ 




1 l^WTS" UTT Stf* Sir 





la^ifl"**" m! .a^ Di atf ab* 44' it* 4b> 



i^ay^4" 5B" aa' *^ 





H'^jy 24* IB* 13* ^ JTO" Sfl* 





1 



1 tf'jS^' 12* 2d* 



aB^aa*'i2" t' sa'^ a*-' 4a" 4z" 



Fig. 11. — i^'-band images of sample galaxies. 
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Offset Category 

Fig. 12. — Plot of the arm/intcrarm ratio in the A'-band image against the offset category 
(1=C, 2=N, 3=A). Points are slightly shifted in the horizontal direction so as not to overlap. 
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Fig. 13. — Plot of A{r, m)/A{r, m = 0) for m = 1 — 4 against the offset category (1=C, 2=N, 
3=A). Points are shghtly shifted in the horizontal direction so as not to overlap. 
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Fig. 14. — Location of the corotation on CO (left) and Ha (right) images. Coordinates are 
(RA, Dec) in J2000. Sohd hne indicates the best-fitted corotation location and dashed lines 
indicate the range of uncertainty. For NGC 0628 and NGC 4254, only dashed lines are drawn 
due to the large radial variation of their rotation curve. 
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Figure [T4l-continued 
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Fig. 15. — Plot of mean H2 surface density S against tgF- Values in the left panel are 
derived by averaging azimuthally, while those in the right are derived by averaging only the 
arm regions. 
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Fig. 16. — Plot of metallicity 12 + log(0/H) versus tsF- Each galaxy has 2 metallicity values 
measured at r = 0.4i?25 (filled circle) and r = 0.8i?s (open circle), where R25 is a radius of 
optical disk and Rg is a scale length. 
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Fig. 17. — Plot of the arm/interarm ratio in the K-hand image against the derived ^sf 
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Fig. 18. — A mean of the relative amplitude A{r,m)/A{r,m — 0) against the derived ^sf- 
Four relative amplitudes of m = 1 — 4 are plotted for each galaxy. Larger values of amplitude 
indicate that the corresponding asymmetric component is larger, presumably due to the spiral 
density waves in that galaxy. 
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